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I . INTRODUCTION 


When  we  model  shock  and  impact  attenuation  we  should  include  nonlinear 
forcing  and  damping  terms1.  After  considering  some  of  the  models  which  have 
been  used  in  the  past,  the  author  studied  some  new  nonlinear  equations  which 
have  applications  in  many  fields  besides  mechanics.  Because  of  their  general 
interest,  some  of  the  solutions  of  one  equation  will  be  described  in  this 
report . 

The  study  of  nonlinear  differential  equations  is  difficult  and  few 
exact  solutions  are  known  for  problems  of  general  interest.  When  exact 
solutions  are  available,  they  offer  several  advantages.  First  of  all,  they 
simplify  tasks  like  mapping  stable  regions,  establishing  the  existence  of 
limit  cycles,  and  so  on.  In  addition,  it  becomes  easier  to  survey  various 
combinations  of  parameters  over  the  entire  ranges  of  their  allowed  values. 
This  can  be  a  fruitful  source  of  ideas.  Finally,  by  avoiding  approximations, 
we  eliminate  an  undesirable  source  of  uncertainty.  With  this  in  mind  we  will 
restrict  ourselves  in  this  report  to  exact  solutions  which  can  be  expressed 
by  a  finite  number  of  elementary  functions. 

The  method  we  will  follow  consists  of  starting  with  an  equation  which 
has  known  solutions,  transforming  it  in  a  general  way  and  then  choosing  a 
particular  form  of  interest.  Since  there  can  be  many  starter  equations, 
transformations  and  particular  forms,  the  case  we  will  discuss  is  merely  an 
example.  The  method  is  not  new,  but  as  we  shall  see,  it  can  lead  to  new 
results  of  general  interest. 

Let  us  begin  with  the  damped,  driven  linear  oscillator 

X  +  M  +  w2X  -  F  (1) 

where  /9  and  u>2  are  constants.  Solutions  and  applications  of  this  equation 
are  well  known.  Next,  we  choose  a  general  transformation  of  the  dependent 
variable  by  considering  X-X[x(t)]  such  that 

X  -  xX'  and  X  -  xX*  +  x2X"  (2) 

where  a  prime  denotes  differentiation  with  respect  to  the  new  variable  x.  Eq 
(1)  becomes 


X'x  +  X"x2  +  £X'x  +  w2X  -  F  (3) 

Finally,  we  choose  a  particular  form 

X-  (  exp (ax)  +  f  (A) 

where  a,  £  and  f  are  constants.  Putting  Eq  (4)  in  Eq  (3)  gives 

x  +  ax2  +  fix  +  (u>2/a)  [  1  +  (f/Oexp(-ax)  ]  -  [ F/(aO  ] exp( -ax)  (5) 
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which  describes  a  forced  nonlinear  oscillator  with  quadratic  as  well  as 
linear  damping.  Because  of  the  forms  in  Eqs  (4)  and  (5)  we  call  this  an 
exponential  oscillator.  Equations  with  linear  rather  than  exponential 
stiffness  have  been  discussed  previously  in  connection  with  impact 
attenuation  using  either  quadratic  damping2  or  combined  Coulomb,  linear  and 
quadratic  damping3.  Equation  (5)  is  not  as  general  as  one  might  like  since 
the  quadratic  damping  coefficient,  a,  appears  in  several  places.  However, 
this  limitation  can  be  advantageous  if  we  desire  a  simple  reduction  to  the 
linear  case.  For  example,  if  £--f-(l/a),  then  x-(l/a)ln(l+aX)  -*  X  and  Eq  (5; 
-*  Eq  (1)  as  a  -»  0. 

Damping  terms  which  depend  on  odd  powers  of  x  like  fix  automatically 
follow  the  sign  of  x.  If  fi>0 ,  this  term  always  opposes  the  motion  (positive 
damping).  If  fi<0 ,  it  always  aids  the  motion  (negative  damping).  However, 
terms  which  depend  on  even  powers  of  x  like  ax2  must  have  the  sense  of  the 
damping  specified.  Various  conventions  have  been  used  for  particular 
applications.  Here  we  will  find  it  more  convenient  to  specify  the  sign  of  o 
for  each  segment  of  a  motion. 

When  the  linear  oscillator  is  externally  excited  by  F(t)  in  Eq  (1) 

we  may  combine  the  right  side  of  Eq  (5)  with  the  last  term  on  the  left  side 

to  make  it  clear  that  the  exponential  oscillator  is  parametrically  or  inter¬ 
nally  excited.  Various  other  forms  of  F  may  also  be  used.  For  example,  if 
F— <pa [ £  exp(ax)  ]-<pa(X-f)  from  Eq  (4),  the  <paX  term  may  be  combined  with  the 
u2X  term  in  Eq  (1),  leaving  ~<pa$  on  the  right  as  a  driving  term.  The  right 

side  of  Eq  (5)  becomes  <p.  If  p  is  a  constant,  the  solutions  are  simple.  If  t p 

is  a  sinusoidal  function  of  time,  Eq  (1)  becomes  a  damped,  driven  Hill's 
equation,  while  Eq  (5)  is  externally  excited,  etc. 

For  constant  F-F0  the  right  side  of  equation  (5)  combined  with  the 
last  part  of  the  left  side  gives  the  stiffness  or  restoring  force 

f  -  (w2/o)(l+6  exp[-Qx])  (6) 

where  6-(f -F0/w2 )/£ .  For  6--1,  a  plot  of  f  versus  x  will  always  pass  through 

the  origin.  This  includes  the  linear  form  w2x  which  results  when  cr-*0.  For 
other  values  of  a,  the  curves  exhibit  "hardening"  or  "softening"  and 
approach  asymptotes  which  depend  on  the  sign  and  magnitude  of  a  for  given 
w2 .  For  6/-1,  the  curves  do  not  pass  through  the  origin.  Different 
parameters  giving  different  curves  may  apply  for  different  segments  of  a 
motion.  For  example,  the  two  curves 


f  -  f(AB)  -  (l-72exp[ -x] ) 
f  -  f (CD)  -  - (l-72exp[x] ) 


(7) 
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with  u2-w2-l,  q-1--q  and  6-6--/2  are  seen  to  be  discontinuous  in  Fig.  1. 

This  restoring  force  will  be  used  later  in  an  example  of  an  oscillation 
between  the  turning  points  X-+.881.  Here  barred  values  apply  for  motion  in 
the  negative  x  direction  while  unbarred  values  denote  motion  in  the  positive 
x  direction.  Another  example  to  be  described  later  uses  the  restoring  force 
segments 


f  -  -1.703(1- .5exp[x] ) 
f  -  1.703(1- .5exp[-x]) 


(8) 


with  u>-1.305,  q-+1  and  5--. 5  which  form  the  closed  loop  ab  in  Fig.  1. 

The  first  f  is  zero  for  x-,693  and  describes  the  lower  segment  while  the 
second  f  vanishes  for  X--.693  and  describes  the  upper  segment.  We  have  not 
used  bars  in  Eq  (8)  to  specify  the  direction  of  motion  for  each 
segment.  A  particular  application  however  will  require  specification  of 
clockwise  or  counterclockwise  motion  around  the  loop.  As  with  any  hysteresis 
loop,  the  integral  of  f  around  the  loop  does  not  vanish. 


Let  us  describe  how  the  loop  in  Fig  1  might  be  used  to  represent  a 

cyclic  stress-strain  ( a-t )  curve  for  example4.  Let  x  represent  a  displace¬ 
ment  which  we  multiply  by  |a|-(l/L)  where  L  is  a  characteristic  length.  Then 
c-|a|x  is  a  strain.  Let  f  represent  a  force  per  unit  mass  which  we  multiply 
by  pL  where  p  is  the  density  of  the  solid  sample.  Then  <7-pLf  is  a  stress. 
Note  that  the  turning  points  of  the  loop  depend  only  on  6  and  a  not  u>.  This 

may  be  seen  from  the  condition  of  loop  closure .namely,  (w2/a) [ l+6exp( -ax0 ) ] 

-- (w2/a) [ l+5exp(ax0) ]  when  we  return  to  the  starting  point  x0 .  This  leads  to 

. 5 [ exp (qx0 )+exp( -ox0 ) ]  -  cosh(ax0)  -  -1/6,  so  if  6--. 5,  then  ax0-+1.317  as 

shown.  Now  suppose  the  sample  we  are  using  is  initially  free  of  stress  and 
strain.  That  is,  we  start  at  the  origin  in  Fig  1  instead  of  at  some  point  on 

the  loop.  For  f-0  when  x-0  in  Eq  (6),  we  must  have  6--1.  Let  pLu7-l  for 
convenience  and  represent  the  first  segment  by  1<r-[ l-exp( - e ) ] ,  using  a-1  and 

«  increasing  (v-x>0) .  In  practice  the  initial  portion  is  a  straight  line  and 
we  should  start  our  first  arc  segment  where  the  elastic  limit  is  reached. 
However,  for  simplicity  we  will  use  an  arc  from  the  beginning  for  purposes 
of  illustration.  In  Fig  2  we  see  that  as  we  increase  the  pull  on  the  sample 
(<7>0  increasing)  we  increase  e  to  1  when  a-. 632  in  the  units  we  are  using. 
Now  suppose  we  begin  to  release  the  tension,  decreasing  e  (v<0).  When  the 
stress  has  been  removed  (a- 0)  we  find  a  residual  elongation  (e>0).  The 
second  arc  segment  is  described  by  2a  -  - [1- . 6exp(«) ]  where  e-|a|x  with  a-- 

1.  The  value  6--. 6  was  found  from  the  continuity  condition  at  the  end  of  the 
first  segment  or  beginning  of  the  second  segment  namely,  7a  -  *[l+6exp(l)]  - 

.632.  To  reduce  the  residual  elongation  to  zero  we  must  apply  pressure 
(o<0) .  If  we  continue  to  increase  the  pressure  after  the  original  length  has 
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Figure  1.  Examples  of  Two-Segment  Restoring  Forces. 


Figure  2.  Approach  to  a  Cyclic  Stress-Strain  Curve 


been  restored,  we  compress  the  sample  («<0),  and  when  2t7--.819,  we  find  e-=- 

1.2.  Now  we  gradually  release  the  pressure,  allowing  the  sample  length  to 
grow  (v>0) .  This  time  when  the  stress  is  removed  we  find  a  residual 
compression,  so  we  must  pull  the  sample  (o> 0)  to  restore  its  original  length 
(e-0) .  Continuity  at  the  beginning  of  the  third  segment  requires  3o 

-( l+6exp [ - ( -1 . 2)  ] )  -  -.819  so  sa  -  [ 1- . 548exp( - £ ) ] .  If  we  pull  until  so  - 

.85  and  e  -  1.3  then  start  to  release,  we  have  4o  -  - [l+6exp(l . 3) ]  -  .851  so 

4a  -  - [ 1- .  504exp ( £ ) ]  and  so  on.  Clearly  S  is  approaching  -.5.  The  final  loop 

is  traversed  repeatedly  and  is  described  by  the  two  segments  a-\ 1- . 5exp( - c ) ] 
and  £7«- [ 1- . 5exp( £ ) ] .  For  a  metal  alloy  sample  the  stress  might  be  measured 

-3 

in  hundreds  of  megapascals  for  microstrains  in  units  of  10  .  Clearly  we  can 

shift  the  position  of  such  loops  as  well  as  change  their  size  and  shape  by 
changing  the  parameters.  Rotations  are  also  possible. 

Consider  another  example  described  by  the  four  segments 

f  -  -8 [ 1 -exp( -2+2x) ] 


f  -  8 [ l-exp(2-2x) ] 

f-  8 [ 1 -exp ( -2-2x) ] 
f  -  -8 [ l-exp(2+2x) ] 


which  have  been  plotted  in  Fig  3  and  might  be  used  to  describe  a  major 
ferromagnetic  hysteresis  loop.  In  the  idealization  of  Fig  3  the  loop  is  not 
quite  closed  for  finite  x.  In  Eq  (9)  w— 4  a— ±2  and  6--exp(+2). 


As  one  may  readily  imagine,  a  motion  may  be  divided  into  more  than  two 
segments,  each  with  its  own  restoring  force,  so  that  a  great  variety  of 
figures  can  be  constructed  by  combining  f-curves.  In  some  cases  we  may 
require  f  and/or  its  derivatives  with  respect  to  x  to  be  continuous  at 
certain  switching  points  in  the  motion  where  parameter  changes  occur.  This 
will  lead  to  relations  between  the  parameters.  If  x  and  t  represent  physical 
variables,  we  may  want  x  and  its  derivatives  with  respect  to  t  to  be 
continuous.  If  X  also  represents  a  physical  variable  instead  of  merely  being 
a  parameter  linking  x  and  t,  we  may  impose  similar  requirements.  In  this 
report  we  will  limit  ourselves  to  motions  with  no  more  than  two  segments  and 
will  usually  consider  x  and  t  (but  not  X)  to  represent  physical  variables. 


If  there  is  no  damping  or  forcing  in  Eq  (1)  it  describes  a 
conservative  system.  In  this  case,  multiplying  by  dX-Xdt-Vdt  and  integrating 
leads  to 

E  -  ,5(V2+w2X2)  -  .  5 (Vq+u?2Xq )  -  . 5 [  (wA  sinwt)2+w*(A  coswt)*]  -  .5(uA)2  (10) 
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Figure  3.  A  Four-Segment  Hysteresis  Curve. 


where  E  is  the  constant  energy  per  unit  mass  in  mechanical  applications. 
Similarly  in  Eq  (5)  with  £-F-0  so  we  can  multiply  by  exp(2ax)dx  with 

dx-xdt  and  integrate  to  find 

2  2  2  2 
e  -  . 5 [ (xexp[ax] )  +  (to/a)  [(exp[ax]+S)  -  S  ]  ] 

-  .  5  [  (x0exp  [ax0  ]  )  +  (u>/a)  [ (exp [ax0 ] +5 )  -  S  ]  ]  (11) 

2 

If  we  multiply  Eq  (11)  by  the  dimensionless  quantity  (q()  we  find 

(a()2e  -  .5[V2  +  (/(X2-f2)]  -  E  -  .5(!»2  (12) 

using  Eq  (4)  and  its  time  derivative  as  well  as  Eq  (10).  Consequently,  there 
can  be  an  arbitrary  constant  energy  in  the  nonlinear  case  even  when  E-0  in 
the  linear  case,  provided  neither  f  nor  w  but  only  A  are  zero.  If  X  changes 
sign  periodically  in  Eq  (4)  then  f  cannot  be  zero  for  real  x.  For  £-1  and 
X-cos(t)  for  example,  we  must  have  f<-l. 

II.  THE  LINEAR  OSCILLATOR 


The  linear  oscillator  is  often  used  to  introduce  many  of  the  phenomena 
exhibited  by  nonlinear  equations,  since,  in  spite  of  obeying  the  principle 
of  superposition,  it  can  also  have  limit  cycle  solutions.  Since  the 
solutions  of  Eq  (5)  are  built  on  those  of  Eq  (1),  it  will  be  convenient  to 
have  some  solutions  of  Eq(l)  exhibited  for  reference. 

It  is  common  practice  to  discuss  solutions  in  the  X,V  phase  plane.  If 
we  let  V-X,  Eq  (1)  becomes 

dV/dX  -  -  (/TV+a>2X- F0 ) /V  (13) 

when  F-F0  is  constant.  Because  of  superposition  we  can  write 

X-X  +  X-  X+  e  (14) 

c  p  c 

where  Xc(t)  is  the  complementary  solution  for  F-0  and  c-F0/w2  is  the 

particular  solution  which  in  this  case  merely  shifts  the  origin.  When 
(0/2) 2>a>2 ,  the  solutions  are  non-oscillatory  with  at  most  one  turning  point 
(V-0) .  The  phase  plane  trajectories  are  given  by 


+ 

[(V-7+Xc)/B+]7  ~  [(V-7'Xc)/B']7 

with 


(15) 
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and 


71  -  -(>3/2)  ±  J(P/ 2)2-u2 


(16) 


B1  -  V0  -  7±Xco  (17) 

for  initial  values  Vn  and  X  —X  -e  by  Eq  (14)  as  shown,  for  example,  by 

co  o 

Minorsky5 .  In  the  Xc , t  plane 

Xc  -  [B'e7+(t'to)  -  B+e7  ]/(7+*7*)  .  (18) 

Special  forms  hold  when  (/?/2 )  2—w2  so  7+— 7  and  B  — B  >  as  Is  well  known. 

When  (£/2)2<w2,  the  solution  is  oscillatory  and  we  have  an  infinite 
number  of  turning  points.  If  we  let  -  (P/2) 2  ,  Eq  (16)  becomes 

7 1  “  -  ( >3/2 )  ±  iwd  (19) 

where  i-/^T  and  u>d  is  the  damped  frequency.  Eq  (18)  with  to-0  becomes 

Xc  -  A  exp(-£/2  t)  cos(u^t-^o)  (20) 

with 

A  -  ±7X2o  +[(V0+^Xco/2)/o;d]2  (21) 

and 

«0  -  arc  tan  [  (V0+£Xco/2)/(W{1Xco) )  .  (22) 

If  we  eliminate  t  between  Eq  (18)  and  its  time  derivative,  we  find  for  the 
trajectory  in  phase  space 

(V+£Xc/2)2  +  (wdXc)2  -  (wdA)2  exp[-<0/«d)(0o-«)]  (23) 

a  logarithmic  spiral.  By  convention  we  may  let  6  Increase  in  a 
counterclockwise  direction  as  t  increases.  The  representative  point  in  phase 
space  then  approaches  or  recedes  from  the  origin,  depending  on  the  signs  of 
P  and  t. 

A  well-known  particular  solution  occurs  when  F  is  sinusoidal.  For 
example,  if  F-F0exp [+iqt ]  where  q-2 nv  is  the  circular  driving  frequency, 
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X  -  (F0/fi) exp  [±iqt  ]  -  (F0/  |fi  | )  [  (cos<f>  +  i  sin^)(cosqt  +  i  sinqt)] 

-  (F0/|fl| ) [cos(qt-4)  ±  1  sin(qt-tf)]  -  Xf  +  i 

-  (F0/ | n | )exp[+i(qt-^) ]  (24) 

with  4>  -  arc  tan  [ q£/(w2 -q2 ) ]  (25) 

and  fi  -  (w2-q2)  +  iq£?  (26) 

with  | Q j  -  700*  where  0  is  the  complex  conjugate  of  0.  If  q  is  zero,  X^  -+ 

F0/w2  as  before.  By  Eq  (14)  the  general  solution  is  the  sum  of  Eqs  (20)  and 

(24).  Since  Eq  (20)  is  a  transient  for  p>0 ,  the  phase  plane  trajectory 
spirals  toward  the  elliptical  limit  cycle 

(qX)2  +  V2  -  (qX0 ) 2  -  2E/m  (27) 

as  X-+Xp  and  the  subscript  p  is  omitted6.  Here  E  is  the  limit  cycle  energy, 
while  m  is  the  mass  or  its  analog  in  a  non-mechanical  application. 

If  F  is  a  square  wave  in  time  instead  of  a  sine  wave  such  that  «-F0/w2 

changes  sign  but  not  magnitude  at  every  turning  point,  we  can  have  another 
example  of  a  limit  cycle5,  provided  £  always  aids  the  motion,  countering  the 
effect  of  p>0 .  Let  us  begin  monitoring  the  motion  at  a  turning  point,  Vo-0. 

Then  by  Eq  (22),  tani90  -  £/(2u>d)  “  c,  defining  c.  By  Eq  (21),  A  - 
X  yi+tan2iS0  -  X  /cosi30  or,  as  in  Eq  (20),  X  -A  cos60-»c,  defining  k. 

CO  CO  CO 

Consequently,  Eq  (14)  becomes 

X0  —  k  +  £  .  (28) 

Since  d— i90-wdt,  then  when  6  has  decreased  by  n  and  0o-d-w  in  Eq 

(23),  while  cos  (wdt^-d0)--cosi90  in  Eq  (20),  then  Xc^--Acosd0  exp(-c»r)  - 
-k  exp(-c?r)  so 

X:  -  £  -  /ce'C*  -  -X0e'C*  +  e(l+e'C")  (29) 

when  we  use  Eq  (28)  to  eliminate  k  in  favor  of  X0 .  This  relates  the  end  of 

the  first  segment  to  the  beginning,  X0 .  For  the  second  segment  (V<0) ,  let  us 

denote  possible  changes  by  using  a  bar  over  every  quantity  except  t.  At  the 
beginning  of  this  segment  t-t^-w/w^,  and  by  analogy  with  Eq  (29) 
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X1  -  e  -  ic  e'C7r  -  -  X0  e’Cir  +  £(l+e'C,r)  (30) 

where  c—p/(2u>^)  .  The  second  form  of  Eq  (30)  comes  from  requiring  X  to  be 

continuous,  namely  X^-X^ .  At  the  next  turning  point  located  on  the  opposite 
side  of  the  origin  from  X^ ,  O  has  decreased  by  2ir ,  so  from  Eq  (20) 

77  —  2  C  7T  ,  «7  ~  ,  “  Cff 

X2-«  +  #ce  -  t  -  (X^-c)e 

-  6  -  [-X0e  +  c(l+e  )  -  e]  e 

-  t ( 1+e  )  -  (  e  (1+e  )  +  X0£  .  (31) 

The  second  form  of  X^  is  obtained  by  eliminating  k  in  favor  of  X^ ,  using  the 

first  form  of  Eq  (30).  The  third  form  of  X2  is  obtained  by  using  the  second 

form  of  Eq  (29),  while  the  final  form  is  simply  a  rearrangement.  If  we  now 
require  X2~£+k  exp[ -2c»r]-X2 ,  X^-e-K  exp  [ -3cir]-X3 ,  and  so  on,  we  can  relate 

each  turning  point  X^  to  X0 .  For  even-numbered  turning  points  we  find 

X2n  “  X°e  [ e ( 1+e  )  -  £  e  (1+e  )]  Z^_0  e  v  '  (32) 

while  for  odd-numbered  points 


X2n+1  -  -e’C"X2n  +  (33) 

where  n-0 , 1 , 2  .  .  .  .  As  the  number  of  turns  increases  without  limit  (n-*®)  ,  the 
first  term  in  Eq  (32)  vanishes  and  the  turning  points  become  independent  of 
the  initial  conditions,  provided  (c+c)>0.  This  illustrates  another  limit 
cycle  solution  of  Eq  (1).  Recall 


„n-l  -k(c+c)*r  r-  -(c+c)*,-! 
Vo  e  -  [1  -  e  v  ] 


Eq  (34)  enables  us  to  write  the  limit  cycle  turning  points,  Eqs  (32)  and 
(33),  in  terms  of  a  finite  number  of  elementary  functions 


X^(even)  -  [«(l+e 


«e*c,r(l+e'c,r)  ]/[l-e 


(c+c)w. 


XL(odd)  -  [  £  (l+e"C,r)  -  £e‘c*(l+e'c,r)]/[l-e'(c+c),r]  .  (36) 
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Eqs  (  35)  and  (36)  could  have  been  found  by  requiring  curve  closure  (X2~X0) 

or  (X^-X^)  as  well  as  continuity.  However,  this  would  not  bring  out  the  fact 

that  the  limit  cycle  is  approached  no  matter  what  the  initial  conditions  may 
be.  In  addition,  continuity  and  closure  also  hold  for  periodic  trajectories 
like  ellipses  which  are  not  limit  cycles. 

If  «-«  and  c-c ,  then  XL(even)-X^(odd)  and  we  have  a  spiral  toward  the 

vertex  e  instead  of  a  limit  cycle.  If  c-c-0  (no  damping)  and  t -t  ,  Eqs  (35) 
and  (36)  are  indeterminate.  The  trajectory  is  an  ellipse  with  center  t  and 
not  a  limit  cycle.  If  c-c-0  and  e/t ,  no  continuity  is  possible.  The 
trajectory  consists  of  discontinuous  elliptical  arcs  and  is  not  a  limit 
cycle.  If  «--£  and  c-c>0 , 


XL(odd)  -  -  XL(even)  -  t  (l+e'C*)/(l-e'C,r)  (37) 

and  the  limit  cycle  is  symmetric  about  the  origin. 

The  restoring  force,  w2X,  in  Eq  (1)  adds  and  subtracts  equal  amounts 
of  energy  each  half-cycle  and  so  is  not  a  factor  in  determining  the 
existence  of  a  limit  cycle,  provided  the  motion  is  oscillatory.  Fig  4  shows 
an  example  ofa  symmetric  limit  cycle  with  - e— €-1/75- . 577 ,  c^w-1.024  and 
/3-/3-.440,  so  w^-w^-1,  c-c-.220  and  0o-arctan  c  -  .217.  From  Eq  (37)  we  find 

X0--73-Xco+ . 577  by  Eq  (14),  so 

X  -  -  ,/S  -  .577  -  -2.309  -  A  cos(-.217)  (38) 

CO 

from  Eq  (20)  and  A--2.364  as  from  Eq  (21).  The  time -dependent  solution  is 
given  by  Eq  (14)  together  with  Eq  (20)  and  its  time  derivative.  For  V>0 

X-  -2.364e' -22tcos(t- .217)  +  .577  (39) 

V  -  2.364e' -22t[ .22cos(t- .217)  +  sin(t-.217)]  (40) 

At  the  first  turning  point  t^-w/w^-ir  and  0^-do-w-  -2.925,  so  cos [ (t-t^) - (d0 - 

*) ]-cos(t- . 217)  once  more.  The  damping  factor  becomes  exp[ - . 22(t-*r) ]  which 
is  equal  to  unity  at  the  beginning  of  the  second  segment  when  t-t^-m.  Then 

X1-75-X  -A  cos ( 2 . 925)  - .  577  ,  so  A-A.  Half-cycles  with  V<0  are  then  described 

by 


X  -  -2.364  e',22(t',r)cos(t-.217)  -  .577 


V  -  2.364  e 


(41) 


- .22(t-*) 


[ ,22cos(t- .217)  +  sin(t- . 217) ] 


(42) 


V 


Figure  U,  A  Linear  Limit  Cycle. 


Examples  of  other  intial  conditions  inside  or  outside  the  limit  cycle 
are  also  shown  in  Fig  4.  If  the  trajectory  starts  from  the  origin  (Xo-Vo-0) , 

we  find  X  -  0  -  .577  -  A  cos(-.217)  so  A  --.591.  Then  for  the  first  half- 
co  o  o 


cycle  with  V>0,  for  0<t<w 

X  -  -.591  e‘ •22(t‘0)cos(t- .217)  +  .577  (43) 

V  -  .591  e' -22(t'0) [ ,22cos(t- .217)  +  sin(t-.217)]  (44) 

Since  X^-X^-. 866-A^cos (2 . 925) - . 577 ,  A^--1.478  and 

X-  -1.478  e' •22(t',r)cos(t-.217)  -  .577  (45) 

V  -  1.478  e‘ •22(t',r)  [  .22  cos(t-.217)  +  sin(t-.217)]  (46) 


for  ir<t<2n .  At  ^2“' 1 • 299-X2-A2Cos ( 6 . 066)+ . 577 ,  so  A2“- 1.921.  If  we  look 

at  the  sequence  -.591,-1.478  and  -1.921,  we  see  that  A  is  approaching  the 
limit  cycle  value , - 2 . 364 ,  even  as  X  and  V  are  approaching  their  limit  cycle 
forms  given  by  Eqs  (39) -(42).  If  the  phase  of  the  driving  force  were  shifted 
initially  by  180°,  the  motion  would  begin  with  X<0,  V<0  and  the  trajectory 
is  a  reflection  through  the  origin  of  the  one  shown.  In  either  case,  the 
system  is  self-excited  and  approaches  a  stable  oscillation  state  even  if 
started  from  rest  at  the  origin. 

If  the  trajectory  starts  from  a  point  outside  the  limit  cycle  such  as 

X0--2,  Vo-0,  we  find  X  --2-.577-A  cos(-.217),  so  A  --2.639  and 
0  ’  0  ’  co  o  o 

X  -  -2.639  e''22(t'0)cos(t-.217)  +  .577  (47) 

V  -  2.639  e' •22(t*°) [ .22  cos(t-.217)  +  sin(t-.217)]  (48) 

for  the  first  half-cycle  with  a  turning  point  at  X^-l . 866-X^-A^cos ( 2 . 925) - 

.577,  so  A^--2.501,  beginning  to  approach  the  limit  cycle  value  of  -2.364. 

The  procedure  is  clear,  so  there  is  no  need  to  elaborate  further  on  this 
example . 


III.  CONSTANT  PARAMETRIC  EXCITATION 


A .  Introduction 

Once  we  specify  0  ,u>7  and  F,  we  can  write  solutions  to  Eq  (1).  If  we 
also  specify  a,£  and  f,  we  can  use  Eq  (4)  to  write  solutions  to  Eq  (5).  The 
number  of  possible  combinations  of  six  constant  parameters  is  too  large  to 
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examine  extensively  in  this  report,  so  we  will  be  highly  selective.  In  this 
section  we  consider  cases  in  which  F-Fq ,  a  constant  (possibly  zero). 

When  F-f-0,  Eq  (5)  is  separable.  Its  solutions  are  well-known  and 
their  application  to  penetration  mechanics  has  been  discussed7. 

B .  Some  non-oscillatorv  trajectories 

Eq  (15)  reduces  to  a  simple  quadratic  form  in  Xc  and  V  when  7+--l  and 

7  --2,  as  occurs  when  w2-2  and  ^-3  in  Eq  (16).  When  X  and  Xc  differ  by  a 

constant  as  in  Eq  (14),  then  the  X,V  trajectory  segments  are  hyperbolic  or 
parabolic  depending  on  the  characteristic  which  in  turn  depends  on  the 

intial  conditions  in  Eq  (17).  Other  values  of  7+  and  7  will  give 
trajectories  which  are  qualitatively  similar. 

In  our  first  example  we  let  w2-2  and  with  X  --1  in  Eq  (17)  so  B+-= 

V0 - 1  and  B  -V0-2.  A  choice  of  V0  in  the  interval  1<V0<2  will  give  B+  and  B 

opposite  signs  so  V,  the  derivative  of  Eq  (18),  will  not  vanish  until  t-**». 

In  other  words,  the  X,V  trajectory  will  have  no  turning  point.  Let  us  take 

V0— 1 . 5  in  appropriate  units.  Then  B+-.5--B  and  the  trajectory  equations  are 

for  q~1 


X  -  -  .Se'^l+e'*)  +  F0/2  -  feX  +  f  (49) 

V  -  e't(e’t  +  .5)  -  f  v  eX  .  (50) 

Clearly  V  and  v  only  vanish  for  t-*»  and  X-*F0/2.  If  Fo<0,  the  X,V  trajectory 

lies  in  the  upper  left  quadrant  of  the  X,V  plane.  If  Fo>0,  it  may  lie 

partially  or  entirely  in  the  upper  right  quadrant.  Other  initial  conditions 
can  lead  to  trajectories  in  the  lower  half  plane.  From  Eq  (50)  we  see  that  v 
has  the  same  or  opposite  sign  as  V  depending  on  the  sign  of  £.  Suppose  £-1 
and  F0-0.  Then  by  Eq  (49),  exp(x0)--l-f>0,  so  f<-l.  Clearly,  the  initial 

value,  x0 ,  as  well  as  the  final  value  can  be  negative,  zero  or  positive.  The 

X,V  and  x,v  trajectories  are  easily  sketched  and  will  not  be  shown  here.  For 
other  values  of  w2  and  f)  the  calculations  are  easier  if  we  work  with 
equations  like  Eqs  (49)  and  (50)  rather  than  with  Eq  (15)  and  its  analog  in 
the  x,v  plane  obtained  by  substituting  Xc«[£exp(ax)+f ] -«  and  V-afvexp(ox)  in 

Eq  (15).  For  applications  to  impact  attenuation,  for  example,  with  given 
initial  conditions  (x0,v0),  we  may  adjust  the  linear  (£) ,  quadratic  (a)  and 
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Coulomb  or  dry  friction  (Fc)  damping  as  well  as  the  other  parameters  in  the 
exponential  stiffness  in  order  to  obtain  a  desired  final  value  of  x  as  v-»0 . 

In  our  second  example,  let  us  consider  a  case  in  which  the  motion  is 
not  only  stopped  but  partially  reversed  (a  rebound  in  a  mechanical  system) . 

If  we  let  X  -0,  then  B+-B  -V  in  Eq  (17).  In  this  case  as  well  as  others  in 
co  o  n 

which  the  V0  term  dominates,  B+  and  B  have  the  same  sign  so  V  as  well  as  v 

will  vanish  after  a  finite  time  and  a  single  turning  point  can  occur  in  each 
of  the  X,V  and  x,v  trajectories.  Here  we  will  consider  some  examples  in 
which  the  turning  points  are  also  switching  points  at  which  the  parameters 
can  change.  In  our  examples  we  will  keep  w2- 2  and  £-3  with  Fo-0  throughout 

the  motion.  If  X0-0,  V0-l  and  -f-f-ar-l,  we  find 


X  -  e*t(l-e’t)  -  ex-l 

(51) 

V  -  e't(2e‘t-l)  -  v  eX 

(52) 

with  the  restoring  force  f-2(l-exp[ -x) ) .  The  first  segment  of  the  X,V 
trajectory  lies  in  the  upper  right  quadrant  as  does  the  first  x,v  segment, 
as  shown  in  Fig  5.  Turning  points  occur  (V-v-0)  when  t^-ln2  at  X^-.25  and 

x^-.223.  If  none  of  the  parameters  change,  then  Eqs  (51)  and  (52)  describe 

the  second  segments  as  well.  Both  of  these  have  minima  followed  by  an 
approach  to  the  origin  as  t-«°.  These  are  shown  by  the  solid  curves  in  the 
lower  half  planes  in  Fig  5.  Different  choices  for  a,£  and  f  could  lead  to 
trajectories  less  similar  than  shown  in  Fig  5.  In  our  example  the  linear 
damping,  the  quadratic  damping  and  the  restoring  force  all  oppose  the  motion 
for  t<t^.  After  t^,  keeping  -f-f-a-1  changes  the  sense  of  the  quadratic 

damping  and  the  restoring  force  from  opposing  to  aiding  the  motion,  while 
the  linear  damping  continues  to  oppose  the  motion. 

Suppose  changes  occur  in  o,(  and  f  at  the  turning  point.  If  x 
represents  a_physical  variable  which  must  be  continuous,  the  condition 
required  is  x^-x^  or 

xx  -  (1/a)  lnKXj-O/?]  -  Xx  -  (1/a)  lnUX^f)/?]  (53) 

where  a  bar  denotes  a  value  after  the  turning  point  and  the  right  subscripts 
denote  values  at  t^.  It  is  possible  to  choose  parameters  which  make  both  x 

and  X  continuous.  However,  if  X  does  not  represent  a  physical  variable,  then 
there  is  no  need  to  do  this.  As  an  example  let  us  take  a--l  so  the  quadratic 
damping  opposes  the  motion  in  the  second  segment  as  well.  The  the  restoring 
force  might  or  might  not  change  its  sense,  depending  on  our  choice  of  f /£ . 
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Figure  5.  Rebound  Trajectorie 


Since  v^-0  at  the  beginning  of  the  second  segment,  f  is  the  only  force 

acting  and  its  sense  influences  the  shape  of  the  second  segment  strongly, 
provided  it  is  not  initially  zero.  If  f^-0,  then  the  x-acceleration  is 

momentarily  zero.  However,  the  motion  continues  because  the  X-acceleration 
is  -  2X^/0  in  our  example. 

Suppose  -f-£-a-  -1  so  f^--2 [ l-exp( . 223) ]-. 5  and  aids  motion  in  the 
negative  x  direction.  Let  us  also  choose  x  continuous  so  X^-X^/(l+X^)- . 2  by 

Eq  (33).  Since  V^-0,  we  find  If*-. 2  and  B  -.4.  The  initial  time  in  Eq  (18)  is 
no  longer  zero  but  t^-ln2.  Then  we  find 


. 8e’t(l-e’t)  -  1  -  e‘x 

(54) 

- .8e't(l-2e*t)  -  v  e‘X 

(55) 

for  the  second  segments.  The  discontinuity  in  X  is  evident  in  the  dashed  X,V 
trajectory  in  Fig  5.  The  new  x,v  trajectory  lies  too  close  to  the  solid  line 
to  be  distinguished  in  the  scale  of  the  figure. 

Suppose  we  choose  f-0.  Then  f--2,  a  constant  opposing  motion  in  the 
negative  x  direction-.  If  we  also  choose  £-.3125,  we  (incidentally)  keep  X  as 
well  as  x  continuous.  The  trajectory  equations  are 

X  -  e ' C ( 1 -e ’ C)  -  .3125  e"X  (56) 

V  -  -e‘t(l-2e't)  -  -.3125  v  e'X  (57) 

so  the  original  solid  line  X^V  trajectory  has  been  restored  as  in  Eqs  (51) 
and  (52).  However,  v>0  when  V<0  by  Eq  (57)  and  the  second  segment  of  the  x,v 
trajectory  is  now  given  by  the  dashed  line  in  Fig  5.  By  dividing  Eq  (57)  by 
Eq  (56)  we  find 

v  -  (l-2e't)/(l-e't)  (58) 

for  t>ln2  and  v-*l  as  t-*»  and  x-*0. 

If  we  include  F0/0 ,  more  possibilities  arise.  For  example,  in  the  last 

case  above,  v>0  can  have  a  maximum  followed  by  an  approach  to  zero  instead 
of  an  approach  to  1  as  in  Eq  (58).  Or  we  can  construct  an  x,v  trajectory 
which  terminates  at  a  particular  value  of  x  after  a  finite  instead  of  an 
infinite  time.  An  example  of  the  latter  situation  occurs  if  we  add  F0--l--F0 
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(a  dry  friction  force  which  always  opposes  the  X  motion).  If  X0-.366  instead 
of  zero  and  V0-l  as  before,  we  find 

X  -  2.732e_t  -  1.866e'2t  -  . 5  -  eX  -  1  (59) 

V  -  -2.732e‘t  +  3.732e*2t  -  v  eX  (60) 

with  as  before  and  x0-.312,  v0-.732  instead  of  0,1.  Now  turning 

points  occur  when  exp ( - t^)-. 732  at  X^-.5  and  x^-,405.  If  a,f  and  £  do  not 
change,  then  X^-X^-.5  by  Eq  (53).  But  F0-l ,  so  X^-X^-,5-0.  Since  V  -0  also, 

we  find  ff^-B  -0  from  Eq  (13).  This  implies  that  X  remains  at  .5  and  V 
remains  zero  with  X-.405,  v-0  for  all  future  time. 

These  examples  should  suffice  to  indicate  some  of  the  types  of  behavior 
which  are  possible. 

C .  Some  Oscillatory  Trajectories 
1 .  Elementary  Examples 

First  of  all,  suppose  0-0,  w2-4  and  F0-0  with  Xo-xo-0  and  V0-v0-l .  In 
addition,  let  -f-{-l/cr~l,  so 


X  -  .5  sin2t  -  eX  -  1 

(61) 

V  —  cos2t  -  v  eX 

(62) 

The  X,V  trajectory  segment  is  an  elliptical  arc,  while  the  x,v  segment  is 
similar  in  appearance.  Both  are  shown  by  solid  lines  in  the  upper  right 
quadrants  of  Fig  6.  The  first  turning  points  occur  at  X^-.5  and  x^-.405  when 

If  we  require  x  to  be  continuous,  we  find  X^-l/3/X^ ,  so  X  is 

discontinuous.  The  second  segments  for  negative  velocities  are  then 
described  by  using  cos [ 2 ( t-x/4) ]-sin2t  so 


X  -  1/3  sin2t  -  1  -  e'X 

(63) 

V  -  2/3  cos2t  -  v  e  X 

(64) 

if  0-0-0,  w-w-2 ,  -f-£-l/a--l.  The  second  turning  points  are  XJ--1/3  and  x^- 
-.288  when  t-3x/4.  To  keep  x  continuous,  we  find  X2~-.25,  so 
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Figure  6.  Examples  of  Oscillatory  Trajectories. 


1 


(65) 


X  -  .25  sin2t  -  eX  - 

V  -  .5  cos2t  -  v  eX  (66) 

describe  the  third  segments,  and  so  on.  It  is  clear  from  Fig  6  that  the  X,V 
trajectory  consists  of  a  series  of  discontinuous  elliptical  arcs  with 
decreasing  amplitudes  1/2,  1/3,  1/4,  etc.,  while  the  x,v  trajectory  is  a 
spiral  approaching  the  origin.  If  a  did  not  change  sign  at  the  turning 
points,  Eqs  (61)  and  (62)  would  describe  the  entire  motion  and  we  would 
obtain  the  closed  trajectories  shown  by  the  dashed  lines  continuing  the 
first  segments.  The  X,V  trajectory  would  be  an  ellipse.  The  alternation  of 
positive  and  negative  quadratic  damping  would  always  bring  the  nonlinear  x,v 
oscillator  back  to  its  original  energy  state.  In  neither  case  do  we  have  a 
limit  cycle. 

If  we  allow  positive  or  negative  linear  damping  we  obtain  spiral  X,V 
arcs  instead  of  elliptical  arcs,  while  the  x,v  trajectory  will  spiral  more 
rapidly  or  more  slowly,  depending  on  the  sense  of  the  linear  damping.  If  we 
also  allow  nonzero  F0 ,  other  possibilities  arise  as  one  can  readily 

imagine . 


2 •  Two-segment  Limit  Cycles 
a .  Introduction 


We  will  limit  ourselves  to  cases  in  which  parameter  changes  occur  only 
at  turning  points  with  a  return  to  the  original  parameter  set  each  time  i9 
changes  by  2m.  We  will  also  keep  x  continuous. 

From  Eqs  (4)  and  (29) 

exp(aXl)  -  C‘1(X1-f)  -  f"1  [  -X0e'C,r+«  (l+e"C,r)  -f  ] 

-  C"1(«-f)(  1+e  C,r)  -  e"C,r  [  ^  *  1(X0  -  C )  ] 

-  b2-e‘C*exp(ax0)  (67) 

where  b2~£  ^(£-f)(l+e  “  -6(l+e  •  Eq  (67)  relates  the  first  turning 

point,  x^,  to  the  initial  turning  point,  x0 .  Continuity  of  x  requires  x^-x^ , 
so  from  Eq  (53) 

Xj_  -  f+?[|'1(X1-f)]Q/Q  -  F+?exp(QXl)  (68) 

using  Eq  (4)  again.  From  the  second  form  of  Eq  (31) 
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(69) 


v  “  C  7T .  n  “  C 

X2  -  t(l+e  )  -  X^e 

Continuity  of  x  requires  *2“x2 ’  so  from  Eq  (53) 

r1(x2-o  -  Kfj^CXj-O]®70  .  (70) 

The  left  side  of  Eq  (70)  is  exp(ax2)  by  Eq  (4).  From  Eqs  (67)  to  (70) 

exp(ax2)  -  [  J~  1(7-jr)  (l+e'C,r)  -  e'c,rexp(ax1)  ]a^a 

r  i  •*  cir  r  x  -cn  /  v  tCi/q  iQ/q  /  -j i  v 

-  [  bx  -  e  [b2-e  exp(ax0)]  ]  (71) 

where  b^-£  ^(c-0(l+e  )  “  -5(l+e  cn) .  Eq  (71)  relates  the  second  turning 

point,  x2 ,  to  the  intial  value,  x0 . 

By  analogy  with  Eq  (67) 

exp(ax^)  -  b2-e  c*exp(ax2) 

-  b2-e'Cir[bre'C,r[b2-e'C%xp(ax0)]Q/Q]Q/Q  (72) 

by  using  the  second  form  of  Eq  (71)  .  By  analogy  with  the  first  form  of  Eq 
(71) 


■  cn 


a/a 


exp(ax^)  -  (b^-e  exp(QX^)] 

-  [bj-e  (b2-e  [b^-e  [h* 


-cn  . .a/a. a/a .a/a, a/a 

-e  exp (ax0 ) ]  )  )  ]  (73) 


using  the  a/a  power  of  Eq  (72),  and  so  on. 

We  can  write  general  expressions  for  even  and  odd  turning  points  if  we 
let  0  and  0  operating  on  a  function  Z  be  defined  by 

OZ  -  b2  -  e'C*Za/a  (74) 

and 

OZ  -  bx  -  e'CnZa/a  .  (75) 

Then  for  the  even-numbered  turning  points 

exp(ax2n)  -  (00)nexp(ax0)  (76) 
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and  for  the  odd-numbered  turning  points 

exp(ax2n+1)  -  Oexp(ax2n)  (77) 

for  n-0,1,2....  The  first  value  of  n  in  Eq  (76)  gives  an  identity,  while  the 
next  two  values  give  the  (a/a)  powers  of  Eqs  (71)  and  (73).  The  first  two 
values  of  n  in  Eq  (77)  give  Eqs  (67)  and  (72).  The  form  of  equations  like 
Eq  (73)  reminds  one  of  a  nest  of  dolls.  Since  a/a  can  have  any  value,  these 
forms  might  also  be  called  generalized  continued  fractions. 

b .  Alternating  Positive  and  Negative  Quadratic  Damping 

In  the  special  case  a/a  -  1,  the  quadratic  damping  is  constant  in 
direction  and  so  alternates  its  sense  from  helping  to  hindering  the  motion 
and  back  every  cycle.  Eqs  (76)  and  (77)  become 

exp(ox2n)  -  e'n(c+c),rexp(QX0)  +  (b^b^)  ^Je‘k(c+c),r  (78) 

and 

exp(ox2n+1)  -  b2  +  aiexp(ax2n)  (79) 

where  we  have  introduced  a^-expC -cir)  and  a2~-exp( -cw) .  Eqs  (78)  and  (79) 

are  the  nonlinear  analogs  of  Eqs  (32)  and  (33).  As  n-*«  both  equations  become 
independent  of  the  initial  condition,  x0 ,  and  we  have  expressions  for  the 

even  and  odd  turning  points  of  a  limit  cycle,  provided  (c+c)>0.  From  Eq  (39) 

exp(axoL)  -  (b1+a2b2)/(l-a1a2)  (80) 

and 

exp(ax1L)  -  (b2+a1b1)/(l-a1a2)  (81) 

where  xqL  is  the  value  of  all  the  even-numbered  turning  points  of  the  limit 

cycle  and  x^  is  the  value  of  all  the  odd-numbered  turning  points.  Eqs  (80) 

and  (81)  are  the  analogs  of  Eqs  (35)  and  (36).  If  the  right  sides  of 
Eqs  (80)  and  (81)  are  negative,  we  need  complex  (ax).  Similar 

considerations  will  arise  in  the  next  subsection  where  we  will  give  an 
example. 

c .  Either  Positive  or  Negative  Quadratic  Damping 

In  some  applications  quadratic  damping  might  have  one  magnitude  and 
sense  (either  always  aiding  or  always  opposing  the  motion).  This  implies  the 
special  case  a/a  —  -1,  so  that  Eqs  (76)  and  (77)  become  continued  fractions. 
For  example,  we  can  write  Eq  (72)  as 
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exp(ax^) -bj  -  a^ 


bl+a2 


exp(ax^) 


(82) 


using  Eq  (67),  relating  the  second  odd  turning  point,  ,  to  the  first  odd 
turning  point,  x^ .  Similarly, 

exp(ax^)-b2  -  a^ 


bl+a2 


b2+al 


b 


,+a 


2 


exp(ox^)  (83) 

and  so  on.  From  Eqs  (82)  and  (83)  we  see  that  we  are  dealing  with  a  periodic 
continued  fraction  with  period  two.  Our  notation  has  been  chosen  to 
facilitate  comparison  with  Wall's  theorem8.  Here  we  will  let  z  represent  the 
limit  of  the  series  of  fractions  begun  in  Eqs  (82)  and  (83).  Wall's 
constants  A^and  are  related  to  our  constants  as  follows:  A^-a^ ,  B^-b^, 

A2~b2a^  and  ®2-b2bl+a2  with  A  ^-B^-l  and  B  ^-Ag-b^-C) .  Then  the  value  of  the 

fraction  (the  limit  of  the  series)  is  one  root  of  the  quadratic  equation 

b^z2  +  (b2b1+a2-a1)z  -  b2a1  -  0  (84) 

In  other  words,  z  -  lim  (exp[ax„  ,i]-b0)  -  exp[ax1T]  -  b„ ,  or 

n-«o  2n+l  2  r  1L  2 


exp[ax1L]  -  [  (b1b2-a2+a1)  ±  7(b1b2+a1-a2) 2  +  4a2bib2  ]A2bi)  (85) 


Since  exp[ax^)  -  bj  +  a^exp[oxQ^]  by  Eq  (67),  we  find 


exp[axQL)  -  t (-b1b2-a2+a1)  ±  7(b1b2+a1-a2) 2  +  4a2bib2^2albl^  (86) 
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Consequently,  the  limit  cycle  odd  and  even  turning  points  represented  by 

and  are  independent  of  the  initial  conditions  and  depend  only  on  the 

parameters.  We  could  obtain  the  same  relations  by  requiring  curve  closure, 
x2-x0,  in  Eq  (71)  or  x^-x^  in  Eq  (72)  or  (82).  Normally,  positive  or 

negative  quadratic  damping  throughout  a  motion  would  lead  to  spirals  toward 
or  away  from  a  fixed  point.  Eqs  (85)  and  (86)  give  the  parameter  relations 
which  are  needed  to  prevent  either  of  these  possibilities  from  happening. 

From  Eqs  (85)  and  (86)  we  see  that  linear  damping  is  not  required  for 
an  x,v  limit  cycle  in  contrast  with  the  requirements  for  an  X,V  limit  cycle. 
Even  if  p-p-0  so  c-c-0  and  a^-aj-'l,  Eqs  (85)  and  (86)  are  still 

determinate.  Positive  or  negative  linear  damping  will  merely  decrease  or 
increase  the  size  of  an  x,v  limit  cycle,  provided  oscillations  are  still 
allowed.  Consequently,  we  can  concentrate  on  examples  which  have  no  linear 
damping  without  much  loss  in  generality.  In  such  a  case  Eqs  (85)  and  (86) 
simplify  to 


exp[ax0L]  -  -6  +  Js2  -  (6/6)  (87) 


exp [ox  ]  -  -6  -  J62  -  (6/6)  .  (88) 

If  we  also  require  symmetry  about  the  origin,  namely  xil“*x0L’  t^ese 

equations  require  6-6<0  with  62>1  for  real  variables  and  parameters.  If 
6<6<0,  the  right  sides  of  Eqs  (87)  and  (88)  will  be  real,  positive  and 
specify  the  turning  points  of  a  real  limit  cycle.  Clearly  there  are  many 
other  possibilities  which  can  make  the  right  sides  of  these  equations 
negative  or  complex  and  require  complex  (ax)  or  complex  x  for  real  a. 

i.  Real  x 


First,  let  us  build  an  x,v  limit  cycle  on  the  X,V  limit  cycle  of  Fig  4 
^although  an  X,V  limit_cycle  is  not  required).  We  recall  that  /9-/9-.44  and 
wc^-1.024  so  u>d-*wd-l ,  c-c-.22  and  0O-.217  with  577,  giving  Eqs  (39)  to 

(42).  Let  us  choose  -a-a-1,  -£-£-1  and  -f— f-»-2.  Then  a2“a^--exp [ - . 22*r]-- . 5 , 

6^—6— -2.577  and  bj-b^-2 . 577(1 . 5)-3 . 866 .  Then  Eqs  (85)  and  (86)  with  positive 

square  roots  give  Xq^--  1.317  and  x^-1.317.  If  we  used  the  negative  square 

roots  we  would  obtain  ±2  which  turn  out  to  be  spurious  (not  the  value  of  the 
continued  fraction).  The  trajectory  equations  become 

X  -  -2 . 364exp ( - . 22t)  cos(t-.217)  +  .577  -  exp(x)  -  2  (89) 

V  -  2 . 364exp( - . 22t) [ . 22cos ( t- ,217)+sin(t-.217)]  -  v  exp(x)  (90) 
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using  Eqs  (39)  and  (40)  together  with  the  nonlinear  forms  of  our  present 
example  for  the  first  half  cycle.  From  Eq  (89)  we  find  that  Xq^--1.317  and 

not  2,  since  Xq^--1.732.  The  segments  with  negative  velocities  are  described 

by 


X-  -2 . 364exp  [-.22(t-ir)  )cos(t  -  .217)  -  .577  -  2  -  exp[-x)  (91) 

V  -  2.364exp[-.22(t-w)][.22cos(t-.217)+sin(t-.217)]  -  v  exp[-x]  (92) 

The  limit  cycles  are  shown  in  Fig  7.  The  X,V  curve  is  the  same  as  the  solid 
curve  in  Fig  4  where  we  illustrated  the  point  that  all  trajectories  approach 
the  limit  cycle  no  matter  what  the  initial  conditions.  Obviously  the  same  is 
true  of  the  x,v  limit  cycle. 


Second,  consider  an  example  without  linear  damping  or  driving  force, 
so  f)-p-e-e-0  and  there  can  be  no  X,V  limit  cycle.  If  -a-a-l-<j-a)  with 


and  -f-f-72,  then  6-6-- Jl.  These  values  in  Eqs  (87)  and  (88)  give  Xq^-.881-= 
*X1L’  us*n6  t*ie  si-Bns  of  the  square  roots  shown.  We  can  reverse  the  signs  of 


the  square  roots  and  still  have  a  single  value  for  the  continued  fraction. 
Therefore,  let  us  use  x^^-- . 881--x^L  instead.  The  restoring  force  is  given 


by  Eq  (7)  and  is  illustrated  in  Fig  1.  As  the  motion  proceeds  from  X--.881 
to  X-.881,  the  force  changes  from  point  A  to  point  B.  At  first  it  is 
negative  and  aids  motion  in  the  positive  x  direction.  Then  it  changes  sign 
and  joins  the  damping  in  opposing  the  motion.  At  the  turning  point  it  is 
discontinuous  and  jumps  to  point  C.  On  the  return  trip  it  first  aids  the 
motion  in  the  negative  x  direction,  then  changes  sign,  and  so  on. 


The  trajectory  segments  with  positive  velocity  are  described  by 


X  -  -cos(t)  -  exp  [x]  -  Jl  (93) 

V  -  sin(t)  -  v  expfx]  (94) 

while  those  with  negative  velocity  are  described  by 

X  -  cos(t)  -  Jl  -  exp[-xj  (95) 

V  -  -sin(t)-  v  exp[-x]  (96) 


The  X,V  trajectory  is  the  unit  circle  while  the  x,v  limit  cycle  is  the 
closed  curve  shown  in  Fig  8.  The  dashed  trajectories  illustrate  how  other 
trajectories  are  attracted  to  the  limit  cycle.  Corresponding  dashed  semi¬ 
circles  are  also  shown. 
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Suppose  -f-f--/26  instead  of  -Jl,  while  all  the  other  parameters 
remain  the  same  as  for  Fig  8.  Then  8—8—-JlZ  and  the  equations  are 


X  - 

-5cos(t)  -  exp[x]  -  Jl 6 

(97) 

V  - 

5sin(t)  -  v  exp[x] 

(98) 

X  - 

5cos(t)  -  JlZ  -  exp[-x] 

(99) 

V  - 

-5sin(t)  -  v  exp[-x] 

(100) 

The  X,V  trajectory  is  a  circle  with  radius  5,  while  the  x,v  limit  cycle  is 

shown  in  Fig  9.  In  Fig  10  we  show  v,t  curves  for  8--Jl  and  -Jl 6.  Such  curves 
were  named  relaxation  oscillations  by  van  der  Pol8. 

Finally,  let  us  consider  a  limit  cycle  which  involves  the  two-segment 
restoring  force  described  by  Eq  (8),  the  loop  in  Fig  1.  As  we  noted  before, 
5-5  must  have  a  value  between  0  and  -1  for  symmetric  f  arcs  concave  toward 
the  origin.  Consequently,  52<1  and  we  cannot  use  Eqs  (87)  and  (88)  for  real 
(ax^)  .  However,  if  /?-/?  is  not  zero  we  can  use  Eqs  (85)  and  (86).  Since 

a^-a2~a--exp ( -C7r)  and  5-5--. 5, we  have  b2~b^-b-. 5(l-a) .Then  c>.838  or  c<-.838 

will  lead  to  real  (ax  ) .  If  c-,838,  (-a)-. 0718,  (b/2)-.268  and  (ax  )-+1.317. 

■L*  L 

If  C--.838,  then  (-a)— 13.91,  (b/2)— 3.728  and  (ax^)- +1 . 317 .  For  either 

equality  the  square  roots  vanish  in  Eqs  (85)  and  (86).  Any  other  value  of  c 
gives  turning  points  different  from  those  shown  in  Fig  1,  which  (as  we 
recall)  were  determined  by  our  choice  of  5  and  a. 

If  we  start  at  xL~  -1.317,  then  v>0  for  the  first  segment  and  a-1 

(positive  quadratic  damping)  means  that  we  follow  the  upper  arc  in  Fig  1 
(the  second  form  in  Eq  (8)  )  and  we  move  clockwise  around  the  loop.  Since 
w2-1.703,  let  us  choose  w^-1  so  c-fi/2-  +.838  and  0o-arctan(c)-+ . 697  for 

V0-0.  In  addition,  let  us  choose  F-0  so  5-f/£--.5.  Then  for  positive  linear 

damping  we  find  X0-A( .  767)--  .  232£  for  x0--1.317  and  at  t-w  we  find  Xj- 

A(  - . 767) ( . 0718)-3 . 232£  for  x^-1.317.  Of  course  there  is  no  A/£  consistent 

with  these  requirements.  However,  if  c-  -.838  (negative  linear  damping), 
then  X0  remains  the  same  while  X^-A( - . 767) (13 . 91)-3 . 232£  and  A/£-  -.303  is 

consistent  with  both  results.  More  generally,  we  find  that  for  quadratic 
damping  of  one  sense  we  must  have  linear  damping  of  the  opposite  sense  in 
order  to  move  clockwise  about  the  loop  of  Fig  1. 
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Figure  9.  Limit  Cycle  for  a  Relaxation  Oscillation. 
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Let  us  choose  -{-{-1  with  the  other  parameters  as  above  so  -f-f«-.5. 
Then  for  a-1  (0<t <n) 

X  -  - . 303exp( . 838t)  cos(t+.697)  -  exp(x)  -  .5  (101) 

V  -  - . 303exp (.838t)[.838cos(t+.697)-sin(t+.697)]  -  v  exp(x)  (102; 
while  for  a--l  (>r<t<2»r) 

X  -  - . 303exp [ . 838(t-w) ]  cos(t+.697  -  -[exp(-x)  -  .5]  (103) 

V-  -  .  303exp[  . 838 ( t-ir)  ]  [ . 838cos ( t+ . 697) -sin( t+. 697) ]  -vexp(-x)  (104) 

In  Fig  11  the  discontinuous  X,V  trajectory  is  shown  by  dashed  curves  while 
the  x,v  limit  cycle  is  shown  by  a  solid  curve.  Our  choice  of  parameters  made 
the  square  roots  vanish  in  Eqs  (85)  and  (86)  and  produced  a  semi- stable  x,v 
limit  cycle  which  attracts  trajectories  from  outside  itself  but  not  from 
inside.  A  sample  attracted  trajectory  is  shown  by  a  dashed  curve  in  rig  11. 
Also  shown  is  a  trajectory  which  begins  at  the  origin  and  soon  becor.-_s 
unstable  (x-«*)  as  do  other  trajectories  which  start  just  inside  the  limit 
cycle  and  begin  to  approach  the  origin. 

Suppose  we  begin  as  before  but  choose  positive  linear  damping  (c-.838) 
and  negative  quadratic  damping  (a--l  implying  the  first  form  of  Eq  (8)  for 
our  first  arc  and  counterclockwise  motion  around  the  loop  in  Fig  1).  Also 
let  F-0  since  this  does  not  affect  the  x,v  cycle.  Then  for  a--l  (0<t<m) 

X  -  -4 . 2 14 exp ( - . 838t)  cos(t-.697)  -  -[exp(-x)  -  .5]  (105) 

V  -  4 . 214exp( - . 838t)  [ . 838cos(t- .697)+sin(t- . 697) ]  -  v  exp(-x)  (106) 

and  for  a-1  (w<t<2m) 

X  -  -4 . 2 14  exp  [  -  .838(t*ir)  ]  cos(t-.697)  -  exp(x)  -  .5  (107) 

V  -  4.214exp[-.838(t-w)]  [ . 838cos(t- . 697)+sin(t- . 697) ]  -  v  exp(x)  (10S) 

The  trajectories  in  Fig  12  are  similar  to  those  in  Fig  11  but  are  canted 
differently.  Now  the  x,v  limit  cycle  attracts  from  the  inside  but  not  from 
the  outside  as  shown  by  sample  dashed  trajectories. 

Of  course  in  applications  like  that  in  Fig  2  there  is  no  limit  cycle. 

ii .  Complex  x 

Let  us  write  x-x^+i  x^,  where  xf  and  x^  are  the  real  and  imaginary 

parts  of  x.  In  addition,  let  us  keep  all  parameters  real.  Eq  (4)  then 
requires  complex  X-Xr+iX^.  For  example,  in  the  simple  case  with  0-F-O,  Eq 

(1)  separates  into  real  and  imaginary  parts  with  solutions 
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ide  Only. 


Figure  12. 


(109) 


X  -  A  cos(wt-0  ) 
r  r  r 

X.  -  A.cos(wt -0. )  (110) 

ill 

Even  with  linear  damping  and  forcing  Eq  (1)  separates  into  independent 
equations  and  both  X^(t)  and  X^(t)  are  easily  written  down.  However,  even 

when  ^-F-0,  Eq  (5)  separates  into  coupled  equations 

x^+  a(x2  -  x2)  +  (w2/a)[l+6  exp( -ax^Jcosax^ ]  —  0  (111) 

x^+  a(2xrx.)  +  (a>2/a)[l-6  exp( -ax^) sinax^ ]  -  0  (112) 

to  these  equations  are  found  by  using  X^_  and  X^  in  Eq  (4)  and 
find 

exp(axr)cosax^  -  (A^/£ ) cos (wt-0^.)  -  8  (113) 

exp(axr)sinax^  -  (A^/f )cos(wt-d^)  (114) 

By  squaring  and  adding  these  equations  we  find  xr(t)  and  by  dividing  one  by 

the  other  we  find  x^(t).  The  same  procedure  may  be  followed  even  when  y9  and 

F  are  not  zero,  although  the  final  expressions  are  more  complicated.  In 
general  we  now  have  four  trajectory  types:  (Xr,Vr),  (X^,V^),  (xr,vr)  and 

(x. ,v. ) . 

l  l 

Consider  the  case  in  which  the  right  sides  of  Eqs  (87)  and  (88)  are 
both  real  but  negative,  for  example  6>0  with  82>8/8.  Since  x  must  be  complex 
these  equations  become 


exP(QXr0L)[coSQXi0L+i 

sinaxiOL]  “ 

-8  +  J87  - 

(8/8) 

(87) 

exp(axrlL) [cosaxilL+i 

sinaxUL]  - 

-s  - 

(6/8) 

(88) 

Since  the  right  sides  are  real,  we  must  have  sinax^-0  at  both  turning 

points.  Since  the  right  sides  are  negative,  we  must  have  cosax, T--l  at  both 

turning  points.  In  other  words,  ox^  must  be  an  odd  multiple  of  k  at  turning 

points.  Since  Eq  (114)  will  then  vanish  at  turning  points,  we  may  choose 
tf^-m/2  if  A^  is  not  zero  and  if  turning  points  occur  whenever  wt  is  a 


The  solutions 
separating  to 


35 


multiple  of  n,  for  example.  If  we  begin  at  a  turning  point  we  can  let  0-0 

in  Eq  (113)  and  compare  with  Eqs  (87)  and  (88)  to  find  (A -  (S/6) . 

In  special  cases  like  {-1/5  or  5-5-1  this  square  root  vanishes,  and 
our  choice  of  O^is  arbitrary. 

Let  0-0  and  0^-w/ 2  in  Eqs  (113)and  (114)  with  w-£-a-l  for 
convenience.  Squaring,  adding  and  taking  the  derivative  gives 


exp(2xr)  -  A?sin2t  +  A2cos2t  -  2Ar5cost  +52  (115) 

v^  -  exp  ( -  2x^_ )  [  (A?  -  A2)cost  +  Ar5]sint  (116) 

while  dividing  and  taking  the  derivative  gives 

tanx^  -  A^sint/tA^cost  -  5)  (117) 

v^  -  A^exp(-2xr)  (Af  -  5cost)  (118) 


when  we  use  Eq  (113)  to  simplify  Eq  (118).  From  the  last  two  equations  we 
see  that  for  A^-0,  v^-0  and  x^  is  the  infinite  set  of  fixed  points  +nn/a 

with  n-0,1,2....  For  non-zero  A^,  5>Ar  allows  an  infinite  set  of  periodic 

solutions,  while  5<Ar  prevents  v^from  vanishing.  From  the  first  two 

equations  we  see  that  there  can  be  only  one  periodic  solution.  Of  course 
unstable  trajectories  of  both  types  are  allowed. 


As  a  particular  example  consider  5-5-/2  with  -A^-A^-JJ*  -  I  -  1  and 


A^-  A^-l .  Then  for  a-1  (0<t<*) 

exp(2xr)  -  3  -  2/2cost 

and  v^  -  72exp(-2xr)  sint 

(119) 

tanx^  -  sint/(cost  -  Jl) 

and  v^  -  exp( -2xr) (l-72cost) 

(120) 

These  equations  describe  the  first  segments  of  the  nonlinear  trajectories. 
Both  linear  trajectories  are  unit  circles  90  degrees  out  of  phase  with  each 
other  in  the  case  we  are  considering.  For  a— 1,  the  second  segments  are 

exp(-2xr)  -  3  +  ijlcost.  and  v^  -  -,/2exp(2xr)  sint 


(121) 


tanx^  -  sint/(-cost  -  Jl)  and  -  exp(2xr)  (-l-72cost)  (122) 

The  (x^v^)  limit  cycle  is  shown  at  the  top  of  Fig  13  while  three  of  the 

(x  v  )  limit  cycles  are  shown  at  the  bottom,  two  as  solid  curves  centered 

at  x^  -  +jt  with  clockwise  motion  of  the  representative  point  as  for  the  real 

trajectory  and  one  as  a  dashed  curve  centered  at  x^-0  with  motion  in  the 

opposite  direction.  The  dashed  trajectory  and  all  those  centered  at  even 
multiples  of  n  are  forbidden  since  cosx^-1  instead  of  -1  as  required  when 

real  turning  points  occur.  The  nonlinear  real  and  imaginary  trajectories  are 
also  out  of  phase.  All  of  these  limit  cycles  are  bistable,  attracting  from 
inside  as  well  as  from  outside.  The  imaginary  trajectories  in  Fig  13 
resemble  the  continents  of  Africa  and  South  America. 

Next,  consider  the  case  in  which  the  right  sides  of  Eqs  (87)  and  (88) 
are  both  complex  because  62<6/6.  By  equating  real  and  imaginary  parts  we 
find 


exp(axr0L)  cosax.0L  -  -6  -  exp(axrlL)  cosax.lL 


(123) 


exp(axr0L)  sinaxiQL  -  JS/1  -  6 2  -  -exp(axrlL>  sinaxilL  (124) 

By  squaring  and  adding  these  equations  we  find 

exp(2axr0L)  -  exp(2axrlL)  -  6/6  (125) 

so  6  and  6  must  have  the  same  sign  and  xr0L-xriL‘  From  (124)  we  also  see 
that  xiOL-XilL'  we  USe  ^  ^^5)  in  Eqs  (123)  and  (124),  we  find 


cosox^  -  ±766  (126) 

for  6<  or  >  0.  By  comparing  Eqs  (113)  and  (123)  we  see  that  for  Ar/0  we  can 
take  19^-19^-71/2  at  turning  points.  By  comparing  Eqs  (114)  and  (124)  we  see 

that  we  can  take  fh-0  with  *  £2- 

Let  6^—ir/2  and  d^-0  in  Eqs  (113)  and  (114)  with  a— 1  again.  If  we 
square  ,  add  and  take  the  derivative  we  find 


37 


(127) 


exp(2xr)  -  A?cos2t  +  A2sin2t  -  2A^5sint  +  52 

vr  -  exp(-2xr) [ (-A?  +  A2)sint  -  A^5  ]  cost  (128) 

while  dividing  and  taking  the  derivative  gives 

tanxi  -  Accost/ (A^sint  -  5)  (129) 

vi  -  Aiexp(-2xr) (-Ar  +  5sint)  (130) 

As  a  particular  example  consider  5-5-1/72  with  A^-A^-0  and  -A^-A^- 
71-.5  -  1/72.  Then  for  ar-1  (0<t<*) 


exp(2xr)  - 

.5(1  +  cos2t) 

and 

V  - 

r 

-.5exp(-2xr)  sint  cost 

(131) 

tanx.  - 
l 

-cost 

and 

vi  - 

.5exp(-2xr)  sint 

(132) 

For  q--1  (n<t<2n) 

exp(-2xr)  - 

. 5(l+cos2t) 

and 

v  - 

r 

-.5exp(2xr)  sint  cost 

(133) 

tanx . - 
l 

-cost 

and 

vi  “ 

.5exp(2xr)  sint 

(134) 

At  the  top  of  Fig  14  the  (xr»vr)  trajectory  is  a  double  egg  with  pausing 

points  (v^-0  without  a  change  in  sign)  at  xro”xri”,°  ■  T^e  turning  points 

xr-±.347  occur  half  way  through  each  segment.  Since  Eq  (125)  is  not 

satisfied,  this  is  not  a  limit  cycle.  The  imaginary  trajectories  at  the 
bottom  of  the  figure  are  not  limit  cycles  either  and  the  middle  one  does  not 

meet  the  condition  cosx^--l/./2  for  5>0. 

Of  course  there  are  many  other  cases  possible  besides  the  simple 
symmetric  examples  we  have  considered.  As  usual,  our  examples  are  meant  to 
illustrate  the  method  rather  than  exhaust  the  possibilities. 


IV.  SINUSOIDAL  PARAMETRIC  EXCITATION 

When  F-F0cosqt  in  Eq  (1),  the  restoring  force  of  Eq  (6)  includes 
sinusoidal  parametric  excitation 
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Figure  14.  Complex 


Periodic  Solutions  without  Attraction. 


f  -  (w2/q) [ 1  +  (l/o [r  -  <cosqt ] exp( -qx)  ] 


(135) 


As  is  well  known,  in  the  absence  of  linear  damping  X(t)  is  the  superposition 
of  free  and  forced  oscillations  and  is  periodic  or  aperiodic  depending  on 
whether  q/o>  is  a  rational  or  irrational  number.  Of  course  none  of  these 
solutions  are  limit  cycles.  When  positive  linear  damping  is  present,  we  have 
the  linear  limit  cycle  of  Eq  (27)  as  the  asymptotic  solution  upon  which  a 
nonlinear  limit  cycle  may  be  built. Since 

x  -  (1/a)  ln[  (X  -  {■)/£]  (136) 

from  Eq  (4),  then  for  real  x,  f  cannot  be  zero  if  X  periodically  assumes 
positive  and  negative  values.  The  stability  of  x  depends  on  the  parameters  f 
and  £  as  well  as  X.  For  £>0,  x  is  stable  provided  f  is  less  than  the  minimum 
value  of  X.  For  £<0,  it  is  stable  if  f  is  greater  than  the  maximum  value  of 
X.  Once  the  parameter  of  X  are  given,  we  can  determine  ranges  of  f  values 
which  lead  to  stable  or  unstable  x.  For  example,  in  Fig  15  we  have  plotted 
maximum  and  minimum  values  of  X  with  zero  damping  and  phase  angle  and 
A-<^-F0-l  so  only  q  varies.  The  curves  shown  connect  all  irrational  values  of 

q/co  and  enable  us  to  obtain  "safe"  ranges  of  f  which  keep  x  stable.  When  q-w 
X  itself  is  unstable.  For  rational  values  of  q/u  the  maxima  and/or  minima  of 
X  may  not  coincide  with  the  curves  shown  as  indicated  by  asterisks  when  q/u 
is  2  or  3.  If  q/u>  is  irrational  and  X  is  aperiodic  the  X,V  trajectory  will 
wander  in  a  limited  region  of  phase  space6.  In  such  a  case  the  x,v 
trajectory  might  appear  to  be  periodic  for  a  time  and  then  suddenly  become 
unstable  when  X-f.'  If  q/u>  is  rational,  X  is  periodic  although  the  trajectory 
might  be  complicated  and  the  period  long.  In  such  a  case  a  stable  x,v 
trajectory  will  also  be  periodic.  However,  if  x(t)  is  observed  for  a  time 
much  less  than  the  period,  it  might  appear  chaotic.  In  practice  it  is 
difficult  to  maintain  frequency  ratios  constant  for  long  times.  In  addition, 
damping  is  inevitable  in  real  systems.  Stability  plots  with  some  other 
parameter  variable  and  the  rest  (including  q)  constant  are  easily 
constructed  as  well. 

Eq  (5)  can  be  written  as  an  autonomous  system  of  three  first  order 
equations.  For  example,  we  may  let  x-v  and  ^>-1  be  two  equations  and 

v  -  -avJ  -  £v  -  (w2/q)[  1  +  (1/0 (  f  +  ecosq^  )exp(-ax)  ]  (137) 

be  our  third  equation.  Some  systems  with  three  degrees  of  freedom,  that  is, 
three  independent  initial  conditions  x0 ,  v0  and  y>0 ,  have  been  found  to  have 

chaotic  as  well  as  periodic  solutions10.  Since  we  have  exact  solutions  for 
the  parametrically  excited  exponential  oscillator,  we  can  distinguish 
periodic  and  aperiodic  solutions  but  find  no  evidence  for  chaos.  If  Eq  (5) 
is  externally  excited,  its  solutions  will  be  based  on  solutions  of  a  Hill's 
or  Mathieu  equation  for  X  and  chaos  may  appear.  Since  these  solutions  cannot 
be  expressed  by  a  finite  number  of  elementary  functions,  they  are  beyond  the 
limits  we  have  set  for  ourselves  here. 
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So  far  in  this  section  we  have  considered  real  forcing  (external  for  X 
and  parametric  for  x) .  If  F  is  complex  then  X  and  x  will  also  be  complex  and 
the  method  of  finding  x(t)  based  on  Eq  (20)  is  similar  to  that  used  in  the 
previous  section.  Instead  of  Eq  (135)  the  restoring  force  will  now  have  a 
term  containing  parametric  forcing  by  a  traveling  wave  exp[i(qt  -  ox.)] .  If 

all  the  parameters  are  real  and  constant  throughout  the  motion,  then  the 
asymptotic  solutions  are 


Xr  -  L  cos (qt-7)  (138) 

X.  -  L  sin(qt-y)  (139) 


where  L“F0/./(u2 -q2 ) 2  +  /32q2  and  tany— /3q/(w2 -q2 )  .  Thus  we  can  find  the 
asymptotic  nonlinear  solutions  from 

£exp(ctxr)  cosctx^  -  Lcos(qt*7)  -  f  (140) 

£exp(axr)  sinax^  -  Lsin(qt-y)  (141) 

From  Eqs  (138)  and  (139)  we  see  that  the  limit  cycle  ellipses  have  the  same 
amplitude  and  frequency  but  are  ninety  degrees  out  of  phase.  Since  any 
parameters  which  allow  stable  oscillations  will  lead  to  limit  cycles,  there 
is  no  need  to  derive  equations  like  Eqs  (85)  and  (86).  From  Eqs  (140)  and 
(141)  we  find 


[£exp(oxr)]2 

-  L2  +  f2  -  2Lfcosr 

(142) 

av 

r 

-  qLf  sinr/[£exp(axr) ] 2 

(143) 

tanax . 

1 

-  L  sinr  /  (Lcosr  -  f) 

(144) 

QVi 

-  qL(L  -  f cosr )  /  [£exp(axr))2 

(145) 

where  r-(qt-7).  If  f  is  zero,  xf  is  a  fixed  point  while  ax^-r-(qt*7)  -»  «  at 
constant  speed.  If  |f|<|L|,  the  (xr.vr)  trajectory  is  a  limit  cycle. 
However,  v^  never  vanishes  so  the  x^  approach  infinity  with  variable  speed. 
If  |  $■  |  —  |  L  |  ,  av^-q/2  independent  of  r  and  all  the  x^  approach  infinity  with 
constant  speed.  The  (xr>vr)  trajectory  is  an  oscillation  between  infinity 
and  a  point  near  the  origin..  If  |f|>|L|,  then  there  are  an  infinite  number 
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of  (x^.v^)  limit  cycles  centered  at  x^  equal  to  multiples  of  it/a  and  one 
(x  , v  )  limit  cycle  which  can  be  symmetric  about  the  origin  if  f2-L2+£2. 


For  example,  suppose  a-L-{-q-l  and  f- Jl .  Eqs  (142)  and  (143)  become 

exp(2xr)-3  -  ijlcosr  and  vr-,/2sinr/exp(2xr) .  Then  for  r  equal  to  even  and 

odd  multiples  of  it  we  have  xrQ--.881  -  -xrl.  Each  imaginary  limit  cycle  also 

has  turning  points  located  symmetrically  on  either  side  of  their  respective 
origins.  However  the  maximum  and  minimum  v^  values  are  not  equal  in 

magnitude  like  the  maximum  and  minimum  v^_  values.  The  trajectories  are 

easily  sketched.  All  are  allowed  and  have  representative  points  which  rotate 
clockwise . 


V.  SUMMARY 

In  this  report  we  have  introduced  the  exponential  oscillator  and 
discussed  some  of  its  exact  solutions.  In  particular  we  have  noted  that 
nonlinear  limit  cycles  can  exist  whether  or  not  linear  limit  cycles  exist, 
at  least  in  some  cases.  Examples  were  given  restricted  to  real  parameters 
but  allowing  complex  variables.  These  were  also  limited  to  one  or  two 
segment  motions.  Included  among  these  examples  was  an  application 
involving  mechanical  hysteresis,  namely  a  cyclic  stress-strain  curve. 

Clearly  we  have  not  exhausted  the  subject  of  the  exponential  oscillator.  We 
merely  mentioned  the  possibility  of  solutions  involving  more  than  two 
segments  as  well  as  solutions  which  cannot  be  expressed  by  a  finite  number 
of  elementary  functions.  It  is  also  clear  that  the  oscillator  we  have 
discussed  is  merely  one  example  of  many  nonlinear  equations  with  exact 
solutions  which  can  be  generated  by  the  same  method.  Not  only  can  we  use 
other  forms  besides  Eq  (4)  in  the  transformation  of  Eq  (3),  we  can  also  use 
other  transformations  of  the  dependent  and/or  independent  variable  singly  or 
in  succession  as  well  as  other  starter  equations.  All  of  these  matters  are 
subjects  for  future  work. 
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